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ABSTRACT 

Digital Images acquired from the geostationary 
METEOSAT satellites are processed and dissemi- 
nated at ESA's European Space Operations 
Centre in Darmstadt, Germany. Their scientific 
value is mainly dependent on their radiometric 
quality and geometric stability. This paper will 
give an overview on the image processing activi- 
ties performed at ESOC, concentrating on the 
geometrical restoration and quality evaluation. 
The performance of the rectification process for 
the various satellites over the past years will be 
presented and the impacts of external events as 
for instance the Pinatubo eruption in 1991 will be 
explained. Special developments both in hard- 
and software, necessary to cope with demanding 
tasks as new image resampling or to correct for 
spacecraft anomalies, are presented as well. The 
rotating lens of MET-5 causing severe geometri- 
cal image distortions is an example for the latter. 


INTRODUCTION 

The history of the METEOSAT system at 
ESA's European Space Operations Centre 
(ESOC) reaches back to the launch of the 
first satellite of the preoperational program in 
1977. While this one was still of an ex- 
perimental nature, its successor became fully 
operational in 1983 and was followed by 


four further satellites since then. MET-3, 
launched in 1988, is the last of the preope- 
rational series and serves to date at 75 °W 
for the Extended Atlantic Data Coverage 
(XADC) mission (de Waard, 1993) on behalf 
of the National Oceanic and Atmospheric Ad- 
ministration (NOAA). MET-4 to MET-6, posi- 
tioned at 0° and 10°W, already belong to 
the operational program (Mason, 1987) and 
are operated on behalf of EUMETSAT. 

An enormous wealth of data has been provi- 
ded by the METEOSAT satellites, with the 
image data being the dominant source of 
information. The images acquired from 
METEOSAT, processed and disseminated at 
ESOC, are used by a wide research com- 
munity as a mean to gain a better under- 
standing of atmospheric processes. This, 
however, depends largely on the radiometric 
quality and geometric stability of the image 
data (Diekmann, 1994). Raw (unprocessed) 
images transmitted from the satellite are dis- 
torted due to its various movements during 
the image taking process. Since this prevents 
an accurate geographical identification of 
image pixels, an evaluation of the image 
distortions and a following resampling of the 
image pixels is a necessary step to reduce 
these geometric errors. This is achieved 
during the onground data handling by an 
image correction process called rectification. 
The principles of this geometrical correction 
scheme is described by Wolff (1985) and 
Bos et al. (1990) summarize the changes 
necessary in order to perform the rectifica- 
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tion process in near real time. This paper 
gives a general overview on the various 
items related to the geometrical processing 
at ESOC. 


METEOSAT IMAGE PROCESSING AT ESOC 


The operational satellites of the METEOSAT 
series are spin-stabilized spacecrafts which 
scan the Earth once every half hour (a slot). 
The detector is a radiometer with a pointing 
direction stepping from south to north by 
rotating the telescope one step per image 
line. Four images are taken simultaneously in 
different frequency bands (infrared window, 
water vapour absorption band, two images in 
the visible spectral region), consisting of 8 
bits per pixel and contain 2500 *2500 pixels 
per image (2500*5000 for each of the VIS 
channels). Two of these data streams are 
currently processed in parallel at ESOC and 
a third one is possible for limited time periods 

(e.g. commissioning, anomaly investigations, 

etc.). A first preprocessing of the continuous 
raw image data stream consists predomin- 
antly of demultiplexing the data to form 
continuous pictures. The bulk of the follow- 
ing image data processing tasks is running 
operational on one of the two mainframe 
(MF) computers (COMPAREX 8/98). The 
second MF serves as backup and is normally 
used for other purposes. Both machines have 
undergone regular upgrades in order to cope 
with the growing need for computing power. 


Figure 1 summarizes the main tasks and 
processes necessary to obtain rectified 
images and meteorological data as the final 
products to be provided to the user com- 
munity. 

After reception and preprocessing of the 
image and corresponding auxiliary data, the 
rectification of the images is the main and 
most time consuming task of the onground 


processing. This term refers to the fact that 
images of the Earth's surface and its atmo- 
sphere transmitted from the satellite are 
distorted due to its various movements 
during the image taking process. Since this 
makes a geographical identification of an 
image element (pixel) almost impossible, an 
evaluation of the image distortion and a 
corresponding resampling of the image pixel 
is a way of reducing these geometric errors. 

The true image signal sensed by the METEO- 
SAT radiometer is additionally altered and 
degraded by various radiometric noise sour- 
ces originating from the satellite and the 
transmission (both Gaussian- and periodic 
noise), by the sampling and digitization 
process, and also through the introduction of 
spatial shift errors (during the attempt to 
geometrically correct the images). A variety 
of parameters is determined during the oper- 
ational image processing for assessing the 
radiometric quality of the METEOSAT images 
(Diekmann and Amans, 1 990; Diekmann 
1994). 

Various monitoring devices serve as indis- 
pensable tools for performance and quality 
controlling at the various stages of the image 
processing chain. Monitors connected to the 
front end processors display different raw 
image channels of the two operational satelli- 
tes before any data processing. Different 
image display and processing devices allow 
an online control of raw and rectified images. 

In addition, a transputer augmented worksta- 
tion (TAW) is available for fast monitoring 
and processing of images. They are trans- 
ferred in real time via an FDDI interface from 
the mainframe computer, reduced to simple 
byte maps and stored on harddisks, which 
are each equipped with a dedicated transpu- 
ter for fast I/O processing. This allows very 
fast and parallel zooming, scrolling and loops 
of up to 500 images in full resolution. A 
semi-automatic software system for quality 
controlling the image rectification (QCIR) is 
run on regular basis. The rectified images are 
finally subject to an advanced segmentation 
process, which is the basis for the 


86 




operational derivation of a number of mete- 
orological parameters in the Meteorological 
Information Extraction Centre (MIEC). 

Special software tools are finally developed 
for satellite commissioning operations and 
special investigations (e.g. satellite anoma- 
lies, end-of-life tests, measurement cam- 
paigns, etc.). 


THE REAL TIME RECTIFICATION PROCESS 

The purpose of the METEOSAT raw image 
rectification is to remove the geometric 
image distortions caused by non-nominal 
spacecraft orbit, attitude, spin and other 
effects. It basically involves the modelling of 
the distortions and the transformation of the 
raw images according to this deformation 
model so as to obtain a rectified image cent- 
red on the nominal sub-satellite point 
(Wolff, 1985). Besides those already men- 
tioned, a number of parameters (radiometer 


stepping parameters, vertical image centre 
positions, interchannel registrations, etc.) are 
used to calculate a pair of deformation matri- 
ces for the horizontal and vertical directions. 

The calculation of some of the deformation 
model parameters is based on a continuous 
analysis of where the southern, northern, 
eastern and western horizons are located in 
the raw IR images. The polar horizon posi- 
tions are necessary for determining paramet- 
ers, such as the vertical image offset, the 
radiometer scanning parameters and for the 
determination of the refined spacecraft at- 
titude. The equatorial horizons, on the other 
hand, are required for the calculation of the 
horizontal centring of the image, the samp- 
ling frequency and for general anomaly 
analysis. "Horizon" in this context is the 
point at which the IR sensor detects a chan- 
ge from space (count < 16) to atmosphere 
(count > 1 6). This IR threshold corresponds 
to a temperature of approximately -90°C. It 
actually reflects a vertical atmosphere co- 
lumn integral and the temperature is the 
mean of this layer at the Earth's limb. Refin- 
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ed polar and equatorial horizons (in fractions 
of a pixel) are finally calculated by fitting a 
second order polynomial to a predefined 
number of lines containing (valid) horizons. A 
warming or cooling of the stratosphere con- 
sequently means a change of the determined 
horizon positions which in return degrades 
the rectification accuracy. This effect is 
observed during the normal seasonal tempe- 
rature changes of the antarctic stratosphere, 
which allowed a modelling of this phenome- 
non, but also after a strong warming due to 
volcanic dust reaching the south polar strato- 
sphere (Diekmann and Bowen, 1992). 

Based on the two deformation matrices, the 
raw image is finally resampled to form the 
rectified image. This process is running since 
several years in near real time (Bos et al., 

1 990), with the old batch-system (starts the 
image rectification after reception of the 
whole image) providing a back-up. In the 
real-time rectification system some of the 
deformation parameters are predicted using 
a combination of information from previous 
IR images and measurements made in the 
first 350 image lines of the current IR image. 

Another prerequisite was the use of the near- 
est neighbour resampling technique, because 
this fast and simple method was the only 
possible for the available data processing 
resources. After many years of experience 
with and improvements of the METEOSAT 
system, the residual geometrical errors in the 
image data are now mainly caused by the 
nearest neighbour resampling (Diekmann and 
de Waard, 1992). Since also the computer 
technologies have vastly stepped forward, 
the use of more sophisticated interpolation 
method within the METEOSAT rectification 
system has become possible. After a tho- 
rough study of various methods, the bicubic 
spline filter was selected for the image re- 
sampling. This CPU time consuming modif- 
ication was first installed and tested on a 
dedicated hardware tool (see below) and 
1993 installed on the operational main- 
frames, whose capacity had been more than 
doubled over the past three years. 
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Figure 2 : Image rectification timing 


The sequence of the main image rectification 
tasks is illustrated in Figure 2. The image is 
received within 25 minutes (5 minutes are 
needed for retrace and a standby period). 
Within the first ca. 4 minutes the deforma- 
tion parameters are determined and the 
deformation matrices calculated. The re- 
sampling starts at that point and catches up 
with the incoming data stream fairly quickly. 
The processed image lines are transferred to 
the dissemination computer and divided into 
dissemination formats. The actual transmis- 
sion of these formats starts with a small 
delay at the end of the slot. In the batch 
processing case this delay is at least eight 
minutes. A real time dissemination system 
with data compression and inclusion of other 
data products in the image data stream was 
developed and tested a few years ago, but 
not used in operations as it necessitates a 
change in the user reception equipment not 
covered by the present programme. 


QUALITY CONTROL OF IMAGE 
RECTIFICATION 

The method used to assess the rectification 
accuracy of METEOSAT images is described 
in detail by Adamson et al. (1988). The 
complex system called "Quality Control of 
Image Rectification" (QCIR) is based on 
about 120 reference landmarks (coastlines, 
islands, lakes) spread over the scanned Earth 
disk. They are extracted from rectified IR and 
VIS images and filtered by a simple automa- 
tic histogram and peak-identification process 



to extract those landmarks with less than ca. 
10% cloudiness. These landmarks are later 
subject of another automatic test (based on 
landmark specific correlation coefficients) to 
delete the remaining cloud contaminated 
landmarks collected during about one week. 
The cloudfree landmarks passing both tests 
are correlated with an accurate digital refer- 
ence landmark data set. The landmark dis- 
placement is defined to be the displacement 
of the maximum of the correlation surface 
from nominal. Results for each landmark are 
presented in terms of line and pixel devia- 
tions as well as in absolute and relative rms 
errors of the sum of both. This process runs 
on a weekly basis. For the relative errors the 
landmark position is compared with the 
results of the previous slot, which gives an 
indication of the rectification stability. 

Constant biases determined with this method 
are usually attributed to a set of registration 


parameters describing the positions of the 
detectors onboard the satellite with respect 
to the detector optical axis. These important 
parameters are usually updated during the 
commissioning of a spacecraft and later 
optimized, if necessary. 

Figure 3 summarizes the performance of the 
METEOSAT image rectification system since 
1989. Large rms errors during the commis- 
sioning of MET-4 in 1989 were caused by 
imperfect detector registration values which 
could be corrected during the following 
weeks. The seasonal wave in the rectifica- 
tion errors can still be identified in 1990; a 
correction scheme based on a model of this 
oscillation was implemented in early 1991, 
resulting in a clear improvement of the qual- 
ity. Volcanic dust in the lower stratosphere 
after the Pinatubo eruption in 1 99 1 caused a 
warming of the east image horizons (July 
and August) and later of the antarctic stra- 
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MET-4 and MET-5 relative landmark correlation results 


rms errors (pixel) 
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Figure 4 : Daily averages of relative rms errors, 1993 - 1994 


better because of the 
higher spatial resolution 
in horizontal direction. 
After start of the bicubic 
spline filter, the geo- 
metrical errors went 
down to values around 
0.2 IR pixels for the 
MET-4 satellite - an im- 
provement of a factor 3. 
A detailed study of the 
impact of the new me- 
thod on the radiometric 
contents proved that no 
series degradations 
lowered the overall 
quality of the rectified 
images. 


tosphere. The resulting rectification degrada- 
tions could be reduced by simple parametri- 
zations of these temperature changes in the 
rectification software {Diekmann and 
Bowen, 1 992). The remaining errors (around 
0.6 IR elements rms) are predominantly in- 
herent in the nearest neighbour resampling. 

In addition to these rather predictable fea- 
tures, the rectification software had regularly 
to be extended in order to cope with unex- 
pected anomalies originating from the space- 
craft. For instance, the "fish" problem of 
MET-4 (electronic disturbances in all image 
channels, Baratelli et al., 1990), and the 
rotating lens of MET-5 (Olivier, 1991) had 
serious impacts on the radiometric and geo- 
metric image quality. 

A significant improvement of the relative 
rectification was achieved in autumn 1993 
with the initiation of the bicubic spline re- 
sampling technique in the operational recti- 
fication system. Figure 4 demonstrates this 
fact. The mean rms errors of landmark dis- 
placements from slot to slot (relative rms) 
was until October 1 993 in the order of 0.6 IR 
pixels - approximately the statistical limit of 
the nearest neighbour technique. The results 
for the VIS channels were always slightly 


THE MET-5 ROTATING LENS ANOMALY 

MET-5 was launched in March 1991 and 
was expected to become the prime operati- 
onal satellite after its successful commis- 
sioning. However, it was discovered that the 
rectification accuracy, in particular the rela- 
tive rectification performance, was signi- 
ficantly worse that the results of other sat- 
ellites. The reason was an unexpected move- 
ment of the Earth's disk in the image frame 
in the order of one IR pixel. Indicator for this 
anomaly were the east/west and north/ 
south horizons. Only the IR and WV channels 
were affected, but not the VIS channels. 

After extensive investigations, the cause of 
the MET-5 image anomaly was identified as 
being caused by a rotating lens (L3) inside 
the radiometer cold assembly just in front of 
the passively cooled longwave detectors 
(Olivier, 1991). This lens is not held firmly 
enough to prevent a rotation. The amplitude 
of this rotation corresponds to geometrical 
distortion of the Earth image of about 1 .1 IR 
pixels, with a frequency between 2 and more 
than 10 slots. Even interruptions of the lens 
rotation have been observed. The more or 
less constant amplitude can be explained if 
one assumes that the optical and geometric 
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centre of the lens do not coincide. Such a 
constant lens rotation introduces a sinusoidal 
distortion into the IR and WV images (be- 
cause L3 focuses the incoming radiation only 
in these detectors) in both horizontal and 
vertical direction. 

As a consequence of this, the image rectif- 
ication system has been extensively modified 
in two steps so as to enable it to minimize 
these geometrical distortions. The software 
version developed first was running in batch 
mode (after reception of the whole image) 
and separates the deformation modelling into 
two parts : the calculation of a base defor- 
mation that corresponds to the deformation 
model that would have been obtained in the 
absence of a lens rotation, and the evalua- 
tion of correction functions to compensate 
for the additional distortions (Hanson and 
Adamson, 1992). The latter is essentially 
achieved by an accurate determination of the 
horizontal Earth horizons and a representa- 
tion of the east-west distortion curve by a 
sine wave model. This distortion curve is 
then translated into the corresponding south- 
north distortion curve using the assumption 
that the lens rotates uniformly during the 
course of a slot. This method gave satisfying 
results, which were, however, still somewhat 
worse than corresponding MET-4 results. 

A real time correction independent on such 
assumptions developed under ESA contract 
is just being installed at ESOC. East and 
west horizons of each line delimit an Earth 
cord. Using datation data which are trans- 
mitted with each image line, the angle bet- 
ween the sun and the midpoint of each chord 
is measured and compared to a predicted 
value. This is based on a METEOSAT state 
vector model determined from observations 
over several days, such that the anomalous 
lens rotation effects are averaged out. When 
the lens motion displaces the line of sight in 
horizontal direction, the measured chord 
midpoint will be shifted with respect to the 
predicted midpoint, which allows a direct 
correction. When the lens rotation displaces 
the line-of-sight in vertical direction, the 
measured angular span of a chord will grow 


or shrink when compared with the predicted 
value based on the state vector. This also 
allows a correction in south-north direction. 
First results have shown that with this me- 
thod a real time rectification is possible with 
a rectification accuracy similar to non dis- 
turbed satellites. 


A TRANSPUTER AUGMENTED 
WORKSTATION 

In spite of various attempts to minimize the 
remaining errors in the real time rectification, 
the nearest neighbour resampling imposed a 
non-reducible limit to the final results, A 
variety of methods were studied which could 
possibly be used instead. All were based on 
interpolation of weighted pixel counts sur- 
rounding the pixel in question. Radiometric 
changes introduced by the finally chosen 
bicubic spline method are well within the 
accuracy of the radiance figures themselves. 
They could even lead to an improvement by 
reducing the impact of the satellite's radio- 
meter and associated equipment on the 
original image. 

The initial technical realization of this very 
time consuming resampling process was part 
of the development of a transputer aug- 
mented workstation (TAW) funded by the 
Austrian Authorities under ESA contract. The 
scope was a computer prototype in hard- and 
software for real time resampling using 
advanced filter techniques, product extrac- 
tion as well as rapid image display and pro- 
cessing. The design of the TAW supports 
interfaces to modules for further applictions. 
Besides the rectification module, a second 
component is connected for near real time 
water vapour wind vectors (WVWV) using an 
optimum pattern matching algorithm. Various 
automatic quality control tools are applied to 
the wind vectors, which were optimized for 
this special application beforehand (de Waard 
et al., 1994). Ad-ditional image processing 
tools are available on a connected worksta- 
tion, which also allows the online calculation 


91 




Figure 5 : TAW functional set-up 


and display of various meteorological pro- 
ducts. 

The TAW is equipped with a total of 23 
CPUs distributed over four different modules 
as sketched in Figure 5 and serves as a 
powerful state-of-the-art development and 
study environment. The capacity of this 
workstation is highlighted by its processing 
power which achieves a performance of 
almost 200 Mflops with the implemented 
applications (Scheiber, 1994). 

The main components of the TAW are the 
resampling and WVWV modules, which both 
consist of T800 transputers and INTEL i860 
processors. These functions are controlled 
via a Sun Sparc host computer, which also 
performs the image data transfer from and to 
the ESOC mainframes via an FDDI interface 
in real time mode. Raw and resampled im- 
ages are also automatically transferred via a 
buffer module to the display system of the 
TAW. Eight harddisks with a total of 3.2 
Gbyte, each equipped with a dedicated 
transputer for fast I/O processing, are avail- 
able to store these images on a cyclic file 
system. This allows fast and parallel zoom- 
ing, scrolling and loops of up to 500 IR im- 


ages on a high resolution monitor. Additional 
functions of this module are pixel inspection, 
window extraction and grey value statistics, 
classifications and Fourier transformations, 
overlays and others. 


SUMMARY 

Digital images in three to four channels of 
two METEOSAT satellites (even three under 
certain circumstances) are received and 
processed at ESOC every half an hour - a 
maximum of 528 images per day. Before 
disseminated to the users, these images are 
geometrically corrected in near real time. The 
quality of this rectification process has con- 
stantly been improved over the past years 
due to better modelling of system inherent 
effects, anomaly handling and hardware 
upgrades. Monthly performance figures for 
perfectly processed and disseminated images 
of in most cases well above 98% and high 
quality cloud motion wind vectors deduced 
from METEOSAT images are good indications 
for a stable and reliable on-ground process- 
ing. Further developments are now coming to 
an end, since EUMETSAT is going to take 
over the whole operations in December 1 995. 
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